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DEFORMATION OF SHALLOW SPHERICAL SANDWICH SHELLS 
UNDER LOCAL LOADING 

By John N. Rossettos 
Langley Research Center 

SUMMARY 

Analytical solutions are obtained for the small symmetrical deformation of shallow 
spherical sandwich caps under both uniform and ring-type local loading with various 
types of edge support. The governing equations, which include transverse shear defor- 
mation, are reduced to a single complex second-order differential equation. The solution 
is obtained in terms of Kelvin functions and a parameter which measures the shear flexi- 
bility of the sandwich core. Explicit formulas for calculating the stresses and displace- 
ments at all points in the cap have been obtained. 

Typical results are presented. They indicate that for localized loading, where the 
radius of the local load is one-tenth of the shell edge radius or less, transverse shear 
deformation plays an important role over a wide range of sandwich core stiffness. In 
general, the stresses and displacements for uniform localized loading are greater than 
those for a corresponding ring load and depend on the loading area. For all cases where 
the edge of the cap is clamped and/or restrained against radial expansion, the maximum 
lateral deflection at the center of the cap varies by at most 20 percent, depending on the 
precise edge condition. The deflections for these cases, however, differ substantially 
(by 60 to 100 percent) from the deflection for the case where the edge is both simply sup- 
ported and unrestrained against radial expansion. 

INTRODUCTION 

The utilization of sandwich-type shells in aerospace vehicles provides structures 
with high strength-weight ratios. One shell component that is widely used is the shallow 
spherical sandwich cap (e.g., the aft heat shield of the Apollo command module and bulk- 
heads for pressurized tankage). Such components are often subjected to various types of 
local loading distributed over a small finite area. Furthermore, the degree of edge 
restraint varies and is largely dictated by the support structure in a specific design con- 
figuration. The present study was undertaken to determine the influence of such factors 
on the spherical sandwich cap. 



References 1 to 4 describe previous studies which involve sandwich-type shells, 
and also include additional bibliography and discussion of past work. Some sample 
numerical calculations, chosen for their simplicity for purposes of evaluating the influ- 
ence of transverse shear deformation in sandwich plates and cylindrical shells, are given 
in reference 2. In general, however, much of the literature applicable to the detailed 
deformation behavior of sandwich shells of revolution has been concerned to a large 
extent with establishing appropriate governing equations for use in applications. (See 
ref. 4.) 

The work herein is concerned with obtaining analytical solutions to the linear equa- 
tions for the symmetrical deformation of spherical caps, taking into account transverse 
shear deformations, as is essential for sandwich shells. These solutions provide explicit 
formulas for calculating the stresses and deformations for all points in the shell. They 
should be directly useful for design as well as for checking large computer programs for 
solution of more complex sandwich-shell problems. An assessment is also made of the 
effects of two types of local load and various edge supports. In order to bracket realistic 
cases of local loading (for example, impact pressure distributions on the Apollo command 
module) two extreme types - a uniformly distributed local load and a circular line (ring- 
type) load — are considered. In conjunction with such loadings, simple and clamped edge 
supports and elastic edge-restraint conditions are treated. 

In the present analysis, the relevant equilibrium equations and boundary conditions 
are derived directly for the shallow spherical shell rather than being obtained as special 
cases of more general theory such as that presented in reference 1. Analytical solutions 
are obtained in terms of Kelvin functions and their first derivatives. These functions 
have been extensively tabulated in the literature. Furthermore, accurate polynomial 
approximations to these functions are available. (See, for instance, ref. 5.) The present 
solutions exhibit certain well-defined singularities as the area of local loading approaches 
zero. These singularities are obtained from the leading terms of the series representa- 
tion of the Kelvin functions for small values of their argument, and reduce to those 
obtained for conventional spherical caps with no transverse shear deformation. 

Calculated values for stresses and displacements which are obtained by means of 
the present solution are given for various conditions of loading and edge support. In 
particular, the effect that the area and type (i.e., uniform or ring) of local loading has on 
deflections and stresses is evaluated. Transverse shear deformation is seen to play an 
important role in the results, especially when the loading is highly localized. This influ- 
ence is depicted by specifying values of a physical parameter which is a measure of the 
transverse shear stiffness of the sandwich shell (this parameter appears explicitly in the 
general solution). The results also indicate important quantitative differences between 
practical types of edge support and provide an understanding of how sensitive the various 
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quantities are to the precise edge condition. Calculations of deflection and moment dis- 
tributions for several areas of localized loading show that only a small portion of the cap 
near the load is sensitive to the precise size of the area of local load. Finally, the role 
of curvature is established quantitatively by obtaining results in a range which includes 
the flat plate as a limiting case. Even slight curvature is found to modify substantially 
the flat-plate values for deflection and bending moment. 

SYMBOLS 



-1 



radius of circle described by cap edge (see fig. 1) 
Kelvin functions 

constant defined by equation (50) 
constant defined by equation (52) 


C, Cj, . . ., c 4 


D ( 


Q 


D 


constants of integration 
radius of circular local load (see fig. 1) 

transverse shear stiffness of isotropic sandwich cap, G c (h + t g ) 


flexural stiffness of isotropic sandwich cap, 


constant of integration 


E s t s h 

<1 - * 2 ) 


d = 


1 + O!^ Ag 


E c 


e s 


sp 


3 W 

DqR 


Young's modulus of face sheets of isotropic sandwich cap 
1 - 1/2 


|l2(l - u 2 }\ 


(see eq. (27b)) 


sandwich core shear stiffness parameter, 


3E„t 


s s 


ft2 / 1111 \j 12(l - i/ 2 )dq a 
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stress function equal to rN r 
complex function, F + Xcp 
shear modulus of sandwich core 
functions defined in appendix B 
horizontal prescribed edge stress resultant 

depth of isotropic sandwich cap measured between middle surface of faces 
functions defined in appendix B 
constants defined by equation (71) 

Kelvin functions 

radial bending moment 
circumferential bending moment 

nondime nsional radial bending moment (see eq. (76)) 

nondimensional circumferential bending moment (see eq. (76)) 

radial stress resultant 

circumferential stress resultant 

nondimensional radial stress resultant (see eq. (76)) 

nondimensional circumferential stress resultant (see eq. (76)) 

magnitude of local load 


intensity of local load in pressure units 



P H 

Py 

Q 

q 

R 

r 

r l’ r 2 

4 s 

u 

V 

V e 

W 

w 

X 



horizontal component of surface loading 
vertical component of surface loading 
transverse shear stress resultant 
nondimensional shear stress resultant (see eq. (76)) 
radius of curvature of spherical cap 
radial coordinate 

values of r at inner and outer edges of a spherical shell 
thickness of faces of sandwich cap 
radial displacement 

vertical stress resultant defined by equation (17a) 
vertical prescribed edge stress resultant 
nondimensional vertical displacement (see eq. (76)) 
vertical displacement 

argument of functions in appendix B, \jjl — 

R 



a 


Y 


€ 


r 



transverse shear parameter, 
transverse shear strain 
radial strain 



- 1/2 
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circumferential strain 


circumferential coordinate 


radial curvature 


circumferential curvature 


Aj, . . Ag functions defined in appendix B 


complex quantity given by equation (27), 


”-;( 3 • e s 2 ) 1/2 M 1 - /2 f 


Poisson's ratio 


rotation of the normal to middle surface of cap 


rotation of the tangent to middle surface of cap 


fij, . . fig functions defined in appendix B 


R i - a 
A S M 1+a 2 


Subscripts: 


annular region, c = r = a (see fig. 1) 


edge value 


e l> e 2 


inner and outer edge of shell (see eq. (7)) 


inner region, r = c (see fig. 1) 


ANALYSIS 


Fundamental Equations 

The governing linear equations for shallow spherical shells undergoing axially 
symmetric displacements, including deformation due to transverse shear, are derived 
with the notation and geometry in figure 1. All quantities (stress resultants, moments, 
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and displacements) are positive as shown. The relevant strain-displacement equations 
are consistent with the Donnell-Mushtari -Vlasov approximation for shallow shells. (See 
ref. 6.) The appropriate strain-displacement relations, which include transverse shear 
deformation (see ref. 1), are 


du r dw 
dr R dr 

(i) 

e = H 
6 r 

(2) 

K = ^ 

r dr 

(3) 

o r 

(4) 

y dr 

(5) 

1! 

(6) 


The rotation of the normal to the middle surface of the shell is denoted by cp, and the 
rotation of the tangent to the middle surface is denoted by \f/. 

The principle of virtual work is written as follows: 



6e r + 6e e + M r 6 * r + M Q 6 k q + Q 6 y )r dr d(9 


p 27T p r 2 p277 

= J J (Py + Ph ^ u ) r d r d# + j l^el + ^el ^ u ) r l ^ 


6w + H 


e2 6u ) r 2 


(7) 


The equilibrium equations and boundary conditions which are a consequence of equa- 
tion (7) and equations (1) to (6) are 


_d_ 

dr 




( 8 ) 
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^( r N r )-N e =-rp H 
A(rM r )-M 0 -rQ = O 


( 9 ) 


( 10 ) 


The boundary conditions at the edges 

1*1 or 

rg are 

to prescribe 


u 

or 

N r 


(11) 

w 

or 

— N 
R r 

-Q 

(12) 

<P 

or 

M r 


(13) 


The stress-strain relations for isotropic sandwich shells (see, for instance, ref. 3) in 
this case become 


«r- A s( N r-'*9) 

(14a) 

£ 9 = A s( N e-^r) 

(14b) 

M r ■ D s( K r + VK n) 

(15a) 

M e - D s( K e * VK r) 

(15b) 

Q = D Q r 

(16) 

Auxiliary quantities, a vertical stress resultant V and a stress function 
defined as follows: 

F, are 

V=iN r -Q 

(17a) 

II 

* 

(17b) 

In the present study, the horizontal component of applied surface force p H 
equal to zero. Equations (9) and (17) then yield 

is taken 

ii 

(18) 
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(19) 


N 


_ dF 
6 dr 



( 20 ) 


The problem can now be reduced to the solution of two coupled second-order dif- 
ferential equations for the variables F and <p. First the vertical equilibrium equa- 
tion (8) is integrated, and equation (17a) is substituted into the result to yield 


rV.-C r Py(r) dr + c (21) 

where c is a constant of integration. From equations (1), (2), and (5) the following 
compatibility equation is obtained: 


e 


r 



( 22 ) 


Equation (22) can be written in terms of F and cp, by substituting from equations (6), 
(14), (16), and (18) to (20), to obtain 


r 2d!F tr dF. F+ rl 
dr 2 dr A s R 


cp - 


R 2 a s d o 


F = - 


r 2 V 

RA s D Q 


(23) 


A second equation relating F and cp is obtained from the moment equation (10) by 
substitution from equations (3), (4), (15), and (20): 


r 2cL| +r d^_ 

dr 2 dr 


<P ~ 


RD, 




(24) 


The two coupled equations (23) and (24) can finally be combined to form a single complex 
second-order differential equation in which the dependent variable is a complex function: 


f = F + \cp 


(25) 


The complex quantity A is chosen so that when equation (24) is multiplied by A and 
added to equation (23) the following equation for f is obtained: 


r^ + r* 
dr 2 dr 


1 +JL M (i + a) 
R 2 


f = -£-p(i + a)V 


(26) 
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The appropriate expression for X is found to be 


where 



( 27 ) 

(27a) 

(27b) 

(27c) 


It is noted that the parameter a (or e g ) is a measure of the influence of transverse 
shear deformation (that is, a — 0 and e g — 0 for no transverse shear deformation). 

General Solution of Governing Equations 

Particular integrals .- A particular integral of equation (26) can be determined once 
the loading Py is specified. For a uniformly distributed load on the cap, Py = p. 

Then, from equation (21), the expression for rV for regions under a uniform load is 

„ 2 

rV = - + c (28) 

2 

With V given by equation (28), a particular solution of equation (26) for regions under a 
uniform load is 


f cR prR 
r “ 2 


(29) 


For regions on which no load acts, p = 0. A particular solution of equation (26) in this 
case is 

f = — (30) 

r 


Homogeneous solution .- The homogeneous solution to equation (26) can be written 
in terms of the modified Bessel functions Ij(\/i+Qf x) and Kj(\/i+a x), in which 

,.*T£ (See ref. 5.) As shown in appendix A, Ij and Kj can be expanded in a 
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series of positive powers of the transverse shear parameter a (values of a are small 
for most practical cases). The coefficients of the various powers of a are expressed 
in terms of Kelvin functions with argument x and their first derivatives. These func- 
tions are extensively tabulated and polynomial approximations of high accuracy have also 
been developed to represent them. (See ref. 5.) The following representations of 
Ij(Vi+o! x) and K^(\/I+a x), as developed and discussed in appendix A, are to be used in 
the present study: 


i^ 2 I^(\/i+a x) = ber' x + -2j-x^ber x - be ^ + ^L_(x 2 bei' x - 2x bei x - 3 ber' x) 


+ 1 


bei’ x+^(beix + be ^ + ^(-x 2 ber’ x + 2x ber x - 3 bei’ x) 


( 31 ) 


-i^^ 2 Kj(^i+Q! x) = ker' x + -^^ker x - — j + ^g-(x 2 kei' x - 2x kei x - 3 ker' x) 

+ i kei’ x + -jp^kei x + + ^g-(-x 2 ker' x + 2x ker x - 3 kei* x) 

(32) 

where primes denote differentiation with respect to the argument x. The homogeneous 
solution to equation (26) is then written as follows: 


f = - ic 2 ) i 1//2 I 1 (^i+o: x) + (c 3 - ic^) -i*/ 2 K^(^i+a x)j 


(33) 


where c^, Cg, Cg, and are real constants. The particular solution (eq. (29)) is 
then added to equation (33) to obtain the complete solution for f . If the definition of f 
(eq. (25)) is taken into account, and real and imaginary parts are equated by using equa- 
tions (29), (31), (32), and (33), the following solutions for the stress function F and the 
rotation of the normal to the middle surface cp are obtained: 


F = c 1 A 1 (q!,x) + c 2 A 2 (Q!,x) + c 3 A 3 (Q!,x) + c 4 A 4 (a,x) 


(34) 


— o — f- <P = c 1 A 5 (a,x) - c 2 A 6 (o!,x) + c 3 A 7 (a,x) - c 4 Ag(o-,x) 
1 + a* 


(35) 
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The functions A 1 (o ,x), . . Ag(a,x) that appear in equations (34) and (35) depend on 
both the parameter a and x, involve Kelvin functions and their first derivatives, and 
are written out in appendix B. 

Stresses and displacements .- The stress resultants N r , N 0 , and Q and the 
moments M r and M 0 can now be obtained by using the expressions for F and cp 
(eqs. (34) and (35)) and previous relations which provide a means of expressing all quan- 
tities in terms of F and cp. The arguments of the functions Aj, Ag, etc., and their 
derivatives with respect to x which appear in what follows are omitted for convenience. 
The results for the stresses are therefore written as 


N 


r 



* 



^4 cR pR 
4 r r 2 2 


(36) 


vr - dF -r '/M a ' , c JE a ' + r ^ A ' + c A*' - - 2E 

N rdF' C lR A l +C 2R A 2 +c 3I A 3 +c 4 R A 4 2 2 


(37) 


Ai Ao Ao A4 

(!s r T= 'lT K !B ,t !R tt 4T 


(38) 


M r = D s(ff +l; f)-( 1+Q!2 ^ ~W S ( c 1^r A 5’ " c 2 V +c 3^r- A 7 ' c 4r A 8 ( 


+ r( C l A 5 " c 2 A 6 + c 3 A 7 " c 4 A 8) 


(39) 


M, 


= D g (l - + j/M r = (l + a^)(l - S ( c 1^5 " c 2 A 6 + C 3 A 7 " C 4 A 8) + 




(40) 


The displacement w can be determined by integrating equation (5). First, if equa- 
tions (6) and (16) are taken into account the expression for can be written as 


dr D, 


Q 


Q 


(41) 
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Equation (41) is then integrated by using equations (35) and (38), and w is written as 


w = (l + a 2 )/iiA s y ( c i A 5 “ c 2 A 6 + C 3 A 7 " 

- j* ( C 1 A 1 + c 2 A 2 + c 3 A 3 + ^4)^ + d 


( 42 ) 


or 


w = (l + a Va s (c 1 H 1 + c 2 H 2 + c 3 H 3 + c 4 H 4) + d 


(43) 


The functions Hj, . . H4 are defined in appendix B, and d is a constant of 
integration. 


SPECIFIC SOLUTIONS 

The general solution is utilized to obtain results for two cases of local loading 
which will bracket other intermediate situations. In one case the local load is uniformly 
distributed over a circular area of radius c with the center at the pole of the spherical 
cap (fig. 1(b)). In the other case a ring load, also of radius c, is applied concentric with 
the cap pole (fig. 1(c)). In each case, various edge-support conditions are treated. 

Uniformly Distributed Local Load 

The boundary conditions which are to be imposed in the case of the uniformly dis- 
tributed local load are now outlined. In order to do this conveniently, two regions of the 
cap are defined. The inner region covered by the local load is defined by r = c, and 
the annular region, where p = 0, is defined by c = r = a. Subscripts i (for inner) 
and a (for annular) will be used to associate each quantity with its respective region. 

In order to avoid repetition in writing equations, it is noted that the expressions for the 
stress function and stress resultants in equations (34), (36), and (37), which hold for the 
inner region, can also be used for the annular region by setting p equal to zero. 

Boundary conditions .- The following conditions are to be applied. At r = 0, 

N r and Ng must be finite, and also 


<Pi = 0 

(44a) 

u i = 0 

(44b) 
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At r = c, the continuity conditions which apply are given by 



V- = V 
v i a 


(45a) 


N r,i = N r,a 


(45b) 


u i = u a 


(45c) 


Wi=w a 


(4 5d) 


<Pi = <P a 


(45e) 

M r,i - M r,a 


(4 5f) 

The boundary conditions at r = a involve the following alternatives: 
Set 1 (no radial restraint, simple support) — 


N =0 
r,a 

M r,a - 0 

o 

II 

a 

* 

(46) 

Set 2 (no radial restraint, clamped) — 




N =0 

r,a 

-s 

p 

1! 

o 

o 

II 

£ 

(47) 

Set 3 (radial restraint, clamped) — 




u a = ° 

o 

II 

w a = ° 

(48) 

Set 4 (radial restraint, simple support) 

— 



O 

11 

3 

M r , a = 0 

o 

II 

oS 

£ 

(49) 


Equations (46) to (49) are actually limiting cases of the more general elastic edge- 
support conditions. In particular, for elastic edge restraint against radial expansion the 
following relation holds at r = a: 



where the constant C u involves knowledge of the edge influence coefficients for the 
supporting structure, which for the present problem can be any shell of revolution. By 
means of equations (2) and (14), equation (50) may be written as 


N 


0,a 



aA s C J Nr ’ a 


= 0 


(51) 


For the case of elastic edge restraint against rotation the following relation holds at 
r = a: 




= a 


(52) 


where the constant also depends on the supporting structure. If equations (3), (4), 

and (15) are used, equation (52) may be written as 


d< ?a 

dr 



i l 

»s 


cp a = 0 


(53) 


It is noted that equations (51) and (53) lead to conditions (46) to (49) when C u and 
take on their limiting values of 0 and °°. 

Evaluation of constants of integ ration .- The conditions at r = 0 are now applied. 
The condition <^(0) = 0 with the use of equation (35) leads to 


lim (cj j. 
x— 0 ' 1,1 


A 5 " c 2,i A 6 + c 3,i A 7 " c 4 J 


,i A s) - 0 


(54) 


where the constants c 1 c 4 . refer to the solution valid in the region r = c. The 

i ,i ,x - 

condition u^(0) = 0, with the use of equations (2), (14), (36), and (37) leads to 


where 


S3, B N e,i “ " N r,i)] = ° 

N fl,i = R^( c l,i A l’ + C 2,i A 2’ + c 3,i A 3’ + c 4,iV ~ ~ ^ 

N r,i = "r ^ c l,i A l + c 2,i A 2 + c 3,i A 3 + c 4,i A 4 + “J - “ 


(55) 

(56) 

(57) 
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It is also required that at r = 0, Ng and N r are to remain bounded. These condi- 
tions, together with equations (54) to (57), which involve limiting procedures as x — 0 
(that is, r — 0), can be evaluated by making use of the expressions for the Kelvin func- 
tions of small argument which are given in appendix B and which appear in the various 
functions needed for the calculation. When this is done, it is found that the following 
constants must vanish: 


c i = 


'3,i 


= c 4,i = 0 


(58) 


The boundary conditions at the edge r = a are considered next, for which the four alter- 
natives are denoted by set 1, set 2, set 3, and set 4. Each set of edge conditions is given 
by three equations which are now discussed. 

Consider the conditions given by set 1 (no radial restraint, simple support). In 
this case equations (46) are to be applied. Use is made of the general solutions 
(eqs. (34) to (43)) with p = 0 since solutions which are valid in the region c = r = a 

are employed here. The argument of the functions involved is J/T (the value of x 

at r = a). For brevity the symbol x & , defined by x & = is employed. Equa- 

tions (46) then lead to the following three equations: 


c a R 


c l,a A l( x a) + c 2,a A 2( x a) + c 3,a A 3( x a) + c 4,a A 4( x a) + — = 0 


(59) 


c l,aS( x a) " c 2,a n 6( x a) + c 3,a°7( x a) “ c 4,a n s( x a) = 0 < 60 ) 

c l,a H l( x a) + c 2,a H 2( x a) + c 3,a H s( x a) + c 4,a H 4( x a) + 5 a = 0 < 61 ) 


The symbols fig, . . ., fig and 

1 d a 

stant d equals ^ — — . 

1+a 2 \/M A s 

to the solution valid in c = r = a. 


Hj, . . ., are defined in appendix B, and the con- 
The constants c i a > • • •> c 4 a > c a > apply of course 


Consider the conditions given by set 2 (no radial restraint, clamped). Two of the 
equations in this case are identical to equations (59) and (61). Equation (60), however, 
must now be replaced by 


c l,a A 5 


( x a) - c 2,a A 6 ( x a) + c 3,a A 7( x a) 


c 4,a A 8( x a 


) = 0 


(62) 
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Consider the conditions given by set 3 (radial restraint, clamped). In this case the 
second and third equations are given by equations (62) and (61). The first equation to be 
used here is 


c l,a n l( x a) + c 2,aS( x a) + c 3,a°3( x a) + c 4,a n 4( x a) " 


1 + V 


c a R ~ 0 


(63) 


Consider the conditions given by set 4 (radial restraint, simple support). The 
three equations to be used in this case are given by equations (63), (60), and (61), 
respectively. 

For each set of the possible edge conditions just described, continuity conditions 
which are given by equations (45a) to (45f) are to be applied at r = c. These conditions 
can be written in terms of the various functions defined in appendix B. First, if equa- 
tion (45a) is applied and use is made of equation (28), the following relation is obtained: 


pc , c i _ c a 
2 c c 

Since c^ = 0, the constant c & can be determined from equation (64) and is 



(64) 


(65) 


where P is the magnitude of the local load. Now, in the five remaining continuity con- 
ditions, equations (45b) to (45f), the argument of the functions which appear is 

For brevity the symbol x-, defined by x- = JJI — , is used in what is to follow. Then, 

c c ’ R 

equations (45b) to (45f) lead, respectively, to the following equations: 


c l,i A l( x c) + c 2,i A 2( x 5) " c l,a A l( x c) ' c 2,a A 2( x 5) - c 3,a A 3( x c) ' c 4,a A 4( x c) = 0 (66 > 

c l,i fi l( x c) + c 2,i°2( x c) " c l,a n l( x c) - c 2,a fi 2( x c) “ c 3,aS( x 5) " c 4,a°4( x c) = ^ 


c l,i H l( x c) + c 2,i H 2( x 5) - c l,a H l( x c) " c 2,a H 2( x c) " c 3,a H 3( x c) " c 4,a H 4( x c) + 5 i " 5 a = 0 

( 68 ) 

c l,i A 5( x 5) " c 2,i A 6 ( x c) “ c l,a A 5( x c) + c 2,a A e( x 5) “ c 3,a A 7( x 5) + c 4,a A s( x c) = 0 ( 69 ) 
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c l,i fi s( x 5) - c 2,i^e( x E) - c l, a H X c) + c 2,a^ 6 ( x c) " c 3,a°7( x 5) + c 4,a fi s( x 5) = 0 


(70) 


The symbols . . ., S2g and Hj, . . . , are defined in appendix B. 

Now, with the constant c & given by equation (65), the remaining eight constants 
of integration, Cj c 2 ^ , c 2 a , Cg a , c 4 a , dj, and d & , can be determined by 

using eight equations. These are the three equations for any set of the edge conditions 
(set 1, set 2, set 3, or set 4) and the five equations ( 66 ) to (70). The values for these 
constants then comprise four different sets of solutions for the case of uniformly distrib- 
uted local loading. The constants are now redefined for convenience as follows: 


c l,i’ c 2,i’ c l,a’ c 2,a’ c 3,a’ c 4,a’ ^i’ ^a “ 77 a ( k l’ k 2’ k 3’ k 4’ k 5’ k 6 ’ k 7’ k 8 ) 


It is observed next that the constants and kg appear in only equations (61) 
and ( 68 ) of the eight equations involved in each set of solutions for the constants. This 
means that only six equations need be solved simultaneously for the six unknown con- 
stants kj, . . ., kg in each case, and the evaluation of only a 6 x 6 determinant is neces- 
sary. Once these six constants are determined, k^ and kg can be calculated by using 
equations (61) and ( 68 ), and they are given by 

- k 8 * k 3 H l(*a) + k 4 H 2( x a) + k 5 H 3 ( x a) + k 6 H 4( x a) < 72 > 

k 7 = k 8 - ( k l - k s) H l( x c) - ( k 2 - k 4) H 2( x c) + k 5 H 3( x c) + k 6 H 4( x c) < 73 > 


For each alternative set of edge conditions the six relevant equations which are to 
be used to determine kj, . . ., kg are as follows: The first two equations for each of 
the sets 1, 2, 3, and 4 are given by equations (59, 60), (59, 62), (63, 62), and (63, 60), 
respectively; the last four equations are the same in each case and are given by equa- 
tions ( 66 ), (67), (69), and (70). 

The matrix equation associated with set 1 (no radial restraint, simple support) for 
calculating kj, . . ., kg is therefore 
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0 

0 

A l( x a) 

A 2( x a) 

A 3( x a) 

A 4( x a) 

k 13 


1 

2 

0 

0 

°5(*a) 

- fi 6( x a) 

°7( x a) 

-°8( x a) 

k 24 


0 

A l( x c) 

A 2( x c) 

0 

0 

- A 3( x c) 

" A 4( x c) 

k 3 


0 

°l( x c) 

°2( x c) 

0 

0 

■°3( x c) 

~ n 4( x c) 

k 4 


a 

c 

A s( x c) 

- A e( x c) 

0 

0 

~ A 7( x c) 

A 8( x c) 

k 5 


0 

fi 5( x c) 

-°6( x c) 

0 

0 

-<M x c) 

°8( x c) 

k 6_ 


0 


where k^g = kj - kg and kg^ = kg ~ k 4> which are used to calculate kj and kg. 

In order to obtain the appropriate matrix equation for the other cases, the following 
changes are to be effected in equation (74). For set 2 (no radial restraint, clamped) the 
six elements of the second row of the square matrix should be replaced by 0, 0, 

Ag(x a ), -Ag(x a ^, A,^x a ), and "Ag^x^, respectively. For set 3 (radial restraint, 
clamped) the first row in the square matrix of equation (74) is to be replaced by 0, 0, 
^l( x a)’ ^2( x a)’ ^3( x a)’ a™* ^4 ( x a) ’ * ke seconc * row is to be replaced by 0, 0, 
Ag( x a ), -Ag^x a ), A^x a ), and -Ag^x^; and the first element of the column matrix on 

the right-hand side of the equation is to be replaced by - i(l + v). For set 4 (radial 

Cl 


restraint, simple support) the first row in the square matrix of equation (74) is to be 
replaced by 0, 0, Og(x a ), ^3( x a)’ anc * ^4( x a)> anc * first element of the 

column matrix on the right-hand side of the equation is to be replaced by - — (1 + v). 

2 


St resses and displacements .- The constants k^, . . ., kg are determined by means 
of equations (72) and (73) and the alternative forms of equation (74) described in the pre- 
ceding paragraph (the constants, of course, take on different values within each case for 
various values of the parameters a^/Rh, c/a, and a). Equation (74) is solved on a 
digital computer in which polynomial approximations of the Kelvin functions, accurate to 
the eighth decimal place (ref. 5), are used. Once the constants are known, the stresses 
and displacements may be determined from equations (34) to (43) by making use of equa- 
tions (58), (65), and (71) and the fact that p is to be set equal to zero in equations (34), 
(36), and (37) when values of quantities are desired in the region c = r = a. In this 
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manner, analytical expressions are determined for the quantities M r , Mg, N r , Ng, 
Q, and w which are valid in each of the two regions r = c and c = r = a. The 
quantities are first nondimensionalized as follows: 


(l - i^ 2 )pR 


n = -IS— N q = — Q 

r ^P r q P W 


27T\//J a 

mr “ (1 + v)PR M r 

a 7ra x T n? 

m e~T » M e W " 


> (75) 


?ra 


(l +a 2 ){JI A PR 


w 

s PR , 


where m r , m^, n r , n^, q, and W are nondimensional. The expressions for m r 


and are then 


m (x) = 


r ' (1 + v)x 


- k 2^6^ 


(r ^ c) (76) 


m r (x) = (l ' + V)x [ k 3^5 (x ) ‘ k 4 fi 6 (x ) + k 5°7 (x ) " k 6°8 (x i] 

Z-TO i p —I 

m e (x) = 2trTv) + 5b A 5 (x) ‘ k 2 A 6 <x U 


(c = r ^ a) 
(r ^ c) 


wn r i n 

m (?( x ) = 2(1 - v) + x C 3 A 5 (X) “ k 4 A 6 (x) + k 5 A 7 (x ) “ k 6 A 8 (x >J ( c = r = a ) 
Expressions for the stress resultants n r , ng, and q are given by 

\2l 


(77) 

(78) 

(79) 


n r (x) = ± 

r X 


n (x) = - 
r v 7 x 


k l A l (x) + k 2 A 2 (x) " 2^(t) 


k 3 Aj(x) + k 4 A 2 (x) + k 5 A 3 (x) + k g A 4 (x) - 


a mr 


n e (x) - k 1 A 1 ’(x) + k 2 A 2 ’(x) - 


(r =1 c) (80) 

(c = r S a) (81) 

(r £ c) (82) 


n 0 (x) = kgAj^x) + k 4 A 2 ’(x) + kgAg'fr) + k g A 4 '(x) 

2Rx 


(c = r = a) (83) 
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q(x) = kjAj(x) + k 2 A 2 (x) 

to 

VII 

(84) 

q(x) = k 3 Aj(x) + k 4 A 2 (x) + k g A 3 (x) + k g A 4 (x) 

(c = r = a) 

(85) 


The vertical displacement W is given by 

W(x) = kjHj(x) + k 2 H 2 (x) + k^ (r § c) (86) 

W(x) = kgH^x) + k 4 H 2 (x) + k 5 H 3 (x) + kgH 4 (x) + k g (5 i r ^ a) (87) 


It is useful to have formulas for computing the vertical displacement W(0) and 
the meridional bending moment m r (0) at the origin, r = 0 (or x = 0). To determine 
the displacement W(0), use is made of equation (86) and of expressions for Hj(x) and 
H 0 (x) for small values of x, which can be obtained from appendix B. It turns out that 

1^(0) = - -y- and H 2 (0) = -1 + so that 

W(0) = k 7 - M kl - (l - ^^)k 2 (88) 


For the bending moment m r (0), use is made of equation (76). In this case the following 
expressions are needed for small x: 


S2 g (x) 

S2g(x) 


i T ±x ( 1 + +°( x2 ) 

> 

x + o(x 2 ) 


(89) 



— — 

1 + V 
\ 2 

k i( x + T')“f k 2 


Equations (89) can be used to evaluate the following limit: 

x-o x " k 2 n 6 (x 3 

If the result in equation (90) is used, equation (76) gives 

m r (0) =^(1 k. 


(90) 


(91) 


which provides a means of calculating the bending moment m r at the origin. 
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Local Ring Loading 

In the case of local ring loading concentric with the pole of the cap (see fig. 1(c)), 
expressions for the stresses and displacements can be obtained by making appropriate 
changes in the previous results for the case of uniformly distributed local load. 

Actually, very few changes are necessary. First, p must be set equal to zero in equa- 
tions (34), (36), and (37) because the vertical component of surface loading is zero every- 
where in the case of a ring load. Another change involves the continuity condition given 
by equation (45a), which must now be replaced by 




(92) 


For the case of ring loading, the conditions at the origin, r = 0, do not differ from those 
for the case of local uniform loading, so that equation (58) applies here also (that is, 
c i = Cg j = C 4 i = 0). Then on using equation (28) with p = 0, equation (92) leads to the 
following result for the constant c & : 


c 


a 


JP 

2rr 


(93) 


The equations which describe the edge boundary conditions, equations (46) to (49) and 
(59) to (63), do not change in the present case of ring loading. Also, the only changes in 
the equations for the matching conditions (eqs. ( 66 ) to (70)) occur in the right-hand sides 
of equations ( 66 ) and (67). After p is set equal to zero in equations (36) and (37), equa- 
tions (45b) and (45c) are used to derive the following two equations, which replace equa- 
tions ( 66 ) and (67), respectively: 


PR 


c l,i A l( x c) + c 2,i A 2( x c) " c l,a A l( x c) " c 2,a A 2( x c) “ c 3,a A 3( x c) " c 4,a A 4( x c) ~ ~ 2rr c 


(94) 


c l,i°l( x c) + c 2,iH X c) - c 1,b.H X c) - c 2,a°2( x 5) " c 3,a fi s( x 5) ' %b.H x c) 


x (1 + y)PR 


27TC 


(95) 


This means that the column vector on the right-hand side of equation (74) is to be 
replaced by 
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1/2 

0 

-a/2c 
(1 + t-)a/2c 
0 
0 


( 96 ) 


and the discussion which follows equation (74) for the various edge conditions applies in 
this case also. Finally, the stresses and displacements can be obtained for the case of 
ring loading by using equations (76) to (87) and by deleting the last term in equations (80) 
to (83). As a result, the formulas for W(0) and m r (0) given by equations (88) and (91) 
also hold for the case of ring loading. In these formulas, the constants kp . . ., kg will 
differ from those for uniform local loading, since the matrix equation used to calculate 
them differs. 


RESULTS AND DISCUSSION 

The analytical solutions derived herein are used to calculate stresses and displace- 
ments for the various cases under consideration. The functions which appear in these 
solutions are written out in appendix B, and the constants of integration are determined 
by evaluating at most a 6 X 6 determinant. This is done for each case of local loading, 
each type of edge support, and several values of core stiffness. 

In this manner, representative results are obtained to show the effects of the type 
of local loading (ring or uniformly distributed load), the area over which it acts, and the 
type of edge support involved. (See figs. 2 to 11.) The results also show the influence 
of the transverse shear flexibility of the sandwich core as measured by values of the 
parameter a (see "Symbols"), which involves the curvature and thickness of the shell 
as well as the transverse shear stiffness of the core. Values for this parameter are 
selected to cover a practical range from strong to relatively weak cores (for honeycomb 
and corrugated cores, for instance, a may range from 0.01 to 0.1). It should be 
pointed out that cases other than those calculated for purposes of the present discussion 
can be studied without difficulty by making use of the solutions within the framework 
developed in this paper. It is noted that when a = 0 the present results reduce to well- 
established results for conventional shells (see, for example, ref. 7). Also, for 


| 

\ 
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isotropic, homogeneous conventional shells with transverse shear deformation, the 
parameter becomes 



Center Deflection and Center and Edge Bending Moments 

The displacement at the center of the spherical cap, w(0), is shown in figures 2 
to 4 for a fixed value of the curvature parameter a^/Rh, which defines the degree of 
shallowness and relative thickness of the shell. In figure 2, w(0) is plotted for various 
values of a against the quantity c/a, which is a measure of the area of local loading. 
The edge conditions are M r = N f = w = 0; that is, simple support and no restraint 
against radial expansion. As the loading becomes highly localized the deflection reaches 
rather large magnitudes, and also the influence of the parameter a becomes more pro- 
nounced. It is noted that for the case of uniform loading over the entire cap surface 
(c/a = 1), values of w(0) for a = 0.001 and a = 0.2 differ by about 20 percent, while 
for c/a = 0.05 the value of w(0) for a = 0.2 is nearly twice that for o; = 0.001. 

Figure 3 shows the effect that different edge conditions have on the deflection at the 
cap center for a given value of a. It is clear that edge restraint is decisive in lowering 
the deflection, since sets 2, 3, and 4, which involve either edge restraint against rotation 
( cp = 0) or edge restraint against expansion (u = 0), have much lower deflections than 
set 1, which involves simple support = 0) and no radial restraint (N r = 0). Figure 4 
shows the influence of the type of loading (ring or uniform) on the deflection. Here 
greater differences in the two curves occur when the load is less localized, as would be 
expected. 

The effects of edge condition, type of loading, and core stiffness on the bending 
moment M r are indicated in figures 5 to 7 for fixed aVRh. In figure 5, M r (0) is 
plotted against c/a for two extreme values of a and two sets of edge conditions for 
uniform loading. The curves indicate that for the larger values of a (i.e., weaker 
core) the bending moment is smaller. This result implies, of course, that the shearing 
action which contributes to the displacement alleviates somewhat the tendency toward 
increased curvature changes during deformation. The effect is small for stiffer sand- 
wich cores. The curves also indicate that for the given boundary conditions, the condi- 
tions which involve greater edge restraint correspond to smaller bending moments. This 
again is due to the relative inhibition of curvature changes. 

Figure 6 demonstrates the effect that the type of loading has on the bending 
moment M r (0). There are substantial differences in the computed values over the 
entire range of c/a, and even for c/a = 0.05 the curves differ by about 20 percent. 
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In figure 7 the bending moment at the edge (r = a) of the shell is given for the fully 
clamped condition in the cases of ring and uniformly distributed loading for various 
values of a. Note that for each value of a, both ring and uniform loading approach the 
same limiting value as the load becomes more localized, while differences occur for 
larger values of c/a. Also, the influence of core stiffness as implied by different values 
of a is shown to be significant for the entire range of c/a. 

Deflection and Bending-Moment Distributions 

In figures 8 and 9 the displacement distribution w(r/a) and bending-moment dis- 
tribution M r (r/a) are plotted for several highly localized loading conditions (c/a = 0.1) 
and several values of a characteristic of typical core stiffnesses. It is noted that in 
figure 8 the influence of the degree of localized loading (values of c/a) on the deflection 
is restricted to a region in the spherical cap which is approximately within a radius of 
r = 0.1a. The effect of different values of the parameter a, however, is felt over a 
greater region of the cap. In figure 9, the major influence of the highly localized loading 
condition on bending moments is again evidenced for r/a = 0.1, with a smaller corre- 
sponding influence of a. 


Effect of Curvature 

In figure 10 the center deflection w(0) caused by a uniform local load over a 
small circular area of radius c = 0.05a is plotted against the curvature param- 
eter a 2 /Rh for various edge conditions. The parameter e cr , is related to a (see 

o / 

"Symbols"), but in contrast to a, it is independent of the curvature parameter a /Rh. 
The solid curves have a value of e ori equal to 0.01, which is representative of honey- 
comb sandwich. For the edge conditions described by sets 1 and 4 (simple support, with- 
out and with restraint against radial expansion, respectively), it can be seen that even 
slight curvature has a significant effect in reducing the flat-plate (a 2 /Rh = o) deflection. 
This is true to a lesser degree for edge conditions given by sets 2 and 3 (clamped, with- 
out and with restraint against radial expansion, respectively). In general, as expected, 
the shells with edges having less restraint experience larger deflections. However, for 
a range of values of a^/Rh between 5 and 10 the deflections w(0) are the same for 
set 3 and set 4. The differences for all cases begin to disappear for shells with high 
curvature. Results corresponding to a weaker core (e g p = 0.02) are given in figure 10 
only for set 1 (although the same conclusions hold for the other edge conditions), and indi- 
cate that the shells of higher curvature are more sensitive to transverse shear effects on 
deflections. The example in figure 10 shows a 7-percent increase in deflection for the 
flat plate (a 2 /Rh = o) and a corresponding 21-percent increase for the deeper shell where 
a 2 / Rh = 20. 
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In figure 11 the meridional bending moment at the edge of the spherical cap (r = a) 
is given for various edge conditions and two values of e . The marked effect of even 

0]J 

slight curvature is again present. It is interesting to note that for the fully restrained 
shell (set 3) the edge moment changes sign as the curvature increases beyond a value 
of a 2 /Rh ~ 6.5. 


CONCLUDING REMARKS 

Analytical solutions are obtained for the linear equations governing the deformation 
of shallow spherical sandwich caps under symmetrical local loading and various edge 
conditions. Results for stresses and displacements typical of a practical range of 
parameter values of edge restraint are presented. Calculations for other cases can be 
made readily, and involve the evaluation of a 6 x 6 determinant. The results indicate 
that for localized loading, where the radius of the local load is 1/10 of the shell edge 
radius or less, transverse shear deformation plays an important role over a wide range 
of sandwich core stiffness. In general, the stresses and displacements for uniform 
localized loading are greater than those for a corresponding ring load. For the cases 
where the edge of the cap is either clamped or restrained against radial expansion, or 
both, the maximum deflection differs by at most 20 percent, depending on the precise 
edge condition. The deflection, however, differs substantially (60 to 100 percent) from 
that for the case where the edge is both simply supported and unrestrained against radial 
expansion. 

The effect of curvature on deflections and bending moments is established by 
obtaining results in a range which includes the flat plate as a limiting case. The results 
show that even slight curvature modifies substantially the flat-plate values for deflection 
and bending moment. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., September 26, 1966, «■ 

124-08-06-17-23. 
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APPENDIX A 


EXPANSION OF BESSEL FUNCTIONS 


In this appendix the modified Bessel functions l^(\[i+a x) and Kj(\/i+ax) are 
expanded in series of positive powers of the small parameter a. The coefficients of the 
various powers of a are expressed in terms of Kelvin functions of argument x and 
their first derivatives. The function 1^ (\/ i+o? x) is treated first, and its series repre- 
sentation as defined in reference 5 (p. 379) leads to 


\2r+l 


r • 2-3 ^ 2r+1 

V 2 / ' 1 rv i m 


h(^)- I-kUrv'- Hr+Tr*#*^ 


r=0 


(Al) 


and 


oo hi ix 




2r+l 


I 1 (\/ix)= V A ?L L = i~-*-/^(ber' x + i bei' x) 

1' A r! r(r + 2) 


r=0 


(A2) 


The term in brackets in equation (Al) can be expanded to give 




2r+l 


1)(- ^ *4 + w) - ^ % 1 . )<2r) 


+ (2r ± l)(2r)( 2 r - _11 io^ 0 ( ffl 4) 

o r G V ' 


(A3) 


Next, the following identities are easily established by using equation (A2): 

- \2r+l 

di^/Tx) = y <2r + ^2 

r=0 


dx 


- (2r + l)(lM 

r! T(r + 2) 


(A4a) 


2 d I^x) ^ (2r + l)(2r)A/Tx\ 
^2 L r! r(r + 2) \ 2 / 
r=0 


r \2r+l 


(A4b) 
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,3 d 3 I 1 (\/Tx) _ y (2r + l)(2r)(2r - 1 ) //Tx \ 2r+1 (A4c 

,3 Lj r! T(r + 2 ) V 2 ) 

** r=0 

When equations (A 3 ) and (A 4 ) are taken into account, equation (Al) leads to the following 
result: 


Ij(\/i+ax) = Ijf/Tx) - 


la , dI l(^~ X ) 
2 dx 




iar 


dljf/ix) 2 d 2 I 1 (\fix) x 3 d^v/Tx) 


16 


dx 


- x 


dx 


dx" 


+ (a 4 terms) 


(A 5 ) 


An expression similar in form to equation (A 5 ) may also be established for the 
function Kj(\[i+oix) by first using its series definition (see ref. 5 ) to write 


K^fi 


l+ax 


)- 4 - 

V IX 


i+h. 


3 a" 

8 


i 5 a" 

16 


0 (a 4 ) 


oo te ^ 2r+1 

y 7 S.B. 

L-i r! 


r=0 


r! r(r + 2 ) r r 


(A6) 


where 


S, = 1 + (2r + 1)(- M + “£ + - (2r + l)(2r)(^ + isgj + (2r + l)(2r)(2r - 1)^ + o(a 4 ) 
B r = log - »[*(r + 1) + *(r + 2)] - if ♦ ^ + M! + 0(a 4 ) 

*(r + 1) = (l + |+|+ • • • +£)- 0.5772157 

and 

» (M 2r+1 

= 7 k + Z r! r/r + 2) B r,0 = -‘' 1/2 (ker' x + 1 tei’ x) (A7) 
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where 


B r,o = log (ir) " l& (r + + + 2 3 


By using equation (A7), the following identities are established in this case: 


dx 


B r,oW^> <A»a) 


-> r=0 


x 


rV^Z) $rV\o * ** - »l(^) < A8b > 

S L J r=0 

Kj (^ix) - — 1-1 = Y (2r ± l)(2r)(2r -.1 + 3x 2 
1 ' VTxJ ^ r! r(r + 2) \2 / r >° ^2 


c 3 

dx 2 


r=0 


dI l('T* x ) / rr \ 

- 3x : - + 21 j (fix) 


dx 


(A8c) 


Equation (A6) can now be expanded by using equations (A7) and (A8). After many terms 
cancel out, the following result is obtained: 


K 1 (\fi+a x) = K x (fTx) - ™ x ^ ^ 


dK^fix) 2 d 2 K 1 (fTx) 


dx 


- x 


dx 


lQf 


16 


x 


dK 1 (^Tx) 2 d 2 K 1 (fTx) x 3 d 3 K 1 (fTx) 


dx 


dx 


3 dx 3 


+ (q:^ terms) (A9) 


Certain relations which can be determined from the properties of Kelvin functions (see 
ref. 5) are to be used next. These are 


ber" x = -bei x - 


ber* x 


bei" x = ber x - 


bei* x 


) (A10) 


ber ,M x = -bei’ x + ber* x bei'" x = ber' x - ^ er -■ + bei' x 

x x 2, x x* 
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ker" x = -kei x 


ker' x 


kei" x = ker x - 


kei* x 


> (All) 


ker"’ x = -kei' x + kei - + \ ker' x kei'" x = ker' x - ker x + kei'x 

X x 2 x x 2 


Finally, equations (A10) and (All) are used together with equations (A2), (A5), (A7), 
and (A9) to obtain 

i 1 / 2 ^ (Ji+a x) = ber* x + y x^ber x - — j + ^g-(x 2 bei' x - 2x bei x - 3 ber ' x) 

+ i bei' x + ^r(bei x + — e - ^’ + ^|-(-x 2 ber' x + 2x ber x - 3 bei’ x) 


(A 12) 


and 


2 

-i 1 / 2 K 1 (\|i+ax) = ker’ x + ^ x^ker x - — x j + ^-(x 2 kei' x - 2x kei x - 3 ker' x) 

+ i kei' x + ^(kei x + - er ' x > ) + — (-x 2 ker’ x + 2x ker x - 3 kei' x) 
2 \ x / 8 


(A13) 


3 

It is noted that the terms containing a are dropped in equations (A12) and (A13) for 
purposes of the present analysis (they may be calculated if desired). For practical 


shallow sandwich shells in the range where 
less than 1 percent is incurred in equations 
of a is truncated to exclude a 2 terms. 


0 < a < 4 
4 


and where 


< an error 
(A12) and (A13) when the series in powers 


a 2 a 2 

Rh 


of 


30 



APPENDIX B 


FUNCTIONS WHICH APPEAR IN THE SOLUTIONS 

The following functions of the parameter a and independent variable x = — 

1 1 

appear in the solutions for stresses and displacements given in the text. They are 
written here for ready reference. 

A 1 (a ,x) = (l + ^)ber' x + (f + ^2?)bei' x+^berx+ ^ bei x 


2 8 / 


4 


A 2 (Q! ,x) = ^1 + ^bei* x - {SL + “^Jber’ x + ^2 bei x - ^2 ber x 
A 3 (a ,x) = (l + ^-)ker’ x + + ^2?Jkei' x+^kerx+ ^2 kei x 


A Aa ,x) = (l + ^)kei' x - k + ^Lr' x + ^ kei x - 


4 

A 5 (a 


ker x 


2 8 jber' x + bei x + ^2ber x 


8 " \2 8 “ 1 2 “ 4 

*V~(§-43 

2\ /„. „ 2..2\ .... 2 


,x> = (i - x + fe . ' " ~ 2 


Ag(a,x) = [l - ^-jber’ x - ^ - ^-JLjbei’ x+^berx- ^j2 be i x 


3a . fai ax \ ker t x + “2 kei x + “12 ker x 


A 7 (a,x)^l--£-jkei’x + (f 


2 8 r 


& 

4 


A 8 (a ,x) = (l - 3|!)ker' ' x - (| - ^)kei’ x + °f ker x - 2^ kei 


AjHa.x) = 
A 2 ’(0! y x) = 
A 3 ’(q!,x) = 


"l 1 " T*) bei x + (“ +£ T") ber x - 1 1 - - (f - 

1 - -^jj-jber x + |o! + — j bei x - ^1 - — 2— + .QL-jkei 2 + ^.2L _ 3o! x jber x 

n< 4 . y _ (t ax ? ■ a^ker* x (a 3o'^x^^kei’ x 

v 8/ \ 2 8/x ” \2 ” 8 / x 


-(1 - — jkei x + 
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Q! 2 \, . / . Q! 2 X 2 


i-T kerx + ri 


keix-fl- — + -)- — + (-- ^-- x2 ) ker ' x 

\ 2 8 / x \2 8 J x 


Iber x + - - bei x - 1 - 


ax 2 3a 2 ^ bei* x (a a 2 x 2 \ber' x 
2 " 8 


/ 2\ 2 2 

-(l-S_)beix + £L5_ 
\ 8 / 8 


Iker x + — x kei x - fl - 


2 

_.„2 


-(l -^Veix+2^ kerx- (l-“2 
\ 8 / 8 \ 2 


3a ^ 

)kei’ x 

8 

1 X 

3a 

2 \ker' : 

8 

1 X 


8 r 

X 

a 2 x 2 ' 

^bei' x 

8 / 

> X 

a 2 x 2 ^ker' x 

8 / 

X 

a 2 x 2 1 

|kei' x 

8 ! 

! 

X 


J 4 (a,x) = ^ Ajdx = ^1 + ^-jber x + + ~ - x j bei x + — bei ' x 

J 2 (a ,x) = J A 2 dx = (l + ^g-jbei x " (^ + ^X“) ber x “ ^ ber ' x 
J 3 (a ,x) = y Agdx = |l + Hp^ker x + + - * j kei x + ^ kei’ x 

J 4 (a, x ) = J A 4 dx =(l + ^-)kei x - (f + ^-jp)ker x “ ^ ker ' X 

> x ) = J A 5 dx = (i - ^-) bei x + (| - ■^x _ ) ber x ber ' x + ^jr bei * x 

J6( a > x ) = j* A 6 dx = - ^j-)ber x ‘ (f " ^T~) hei x + ^Y bei ' x +9L y ter' x 

J 7 (a,x) = J Arydx = ^1 - ^-jkei x + ^ - ~ g X - j ker x - ^ ker’ x + kei' x 

J 8 (a ,x) = J Agdx = (l - ^-jker x " (f " '^ L_ ) kei X+ ^T kei ' x + ker ' X 
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Definitions : 


= xAj T - t'Aj 

S = xA 5* + ^5 

= xAg* " v ^2 

= xA0 + J/A g 

fig = XAg' - tAg 

Q,rj = XA 7 ? + VKrj 

4^ 

1! 

£ 

4^ 

£ 

4^ 

fig — xAg + tAg 

Hj = Jg - 2aJ ^ 

Hg — Jiy “ 2(jJg 

h 2 = ~ J 6 - 2q;J 2 

H4 = -Js " 2 qjJ4 


For small values of the argument of the functions ber, bei, ker, and kei, the following 
series are useful (ref. 5): 


ber x = 1 


2 2 • 4 2 


+ . • . 


beiX= V“T^2 2 

q 2 Z . 4 ^ • 6 ^ 


+ . . . 


ber’ x = - 5- + . . . 
16 


bei' x = ^ - 


X 


2 2 2 • 4 2 • 6 


ker x = -In x + 0.1159 + + . . . 

2 2 

kei x = -2— lnx-^ + 1.1159 

4 4 4 


ker' x = - — + ^2£ + 2_lnx + . . . 
x 8 16 


kei' x = -- lnx-- + 0.558x + . . . 
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c/a 

Figure 6.- Comparison of bending moment at center of cap for uniform and local ring loading. a'v/Rh = 2.5; a = 0.2; v = 0.3. 




Figure 7.- Meridional bending moment at edge of cap for uniform and local ring loading where the edge is clamped. a 2 /Rh = 2.5; 




Figure 8.- Displacement w(r/a) caused by load P uniformly distributed over small circular areas defined by three values of c/a. 



Figure 9.- Meridional bending moment caused by load P uniformly di: 



r/a 


over small circular areas defined by three values of c/a. 




•w(O) 



a 2 /Rh 


Figure 10.- Center deflection due to uniform local load (c/a = 0.05) for various edge conditions plotted against curvature parameter a 2 /Rh. 

v = 0.3. 
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